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Abstract 

In this paper we investigate a stochastic model for an economic 
game. To describe this model we have used a Wiener process, as the 
noise has a stabilization effect. The dynamics are studied in terms 
of stochastic stability in the stationary state, by constructing the 
Lyapunov exponent, depending on the parameters that describe the 
model. Also, the Lyapunov function is determined in order to ana- 
lyze the mean square stability. The numerical simulation that we did 
justifies the theoretical results. 
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1 Introduction. 



Stochastic modeling plays an important role in many branches of sci- 
ence. In many practical situations, perturbations are expressed in terms of 
white noise, modeled by brownian motion. The behavior of a deterministic 
dynamical system which is disturbed by noise may be modeled by a stochastic 
differential equation (SDE), [8j. The stochastic stability has been introduced 
by Bertram and Sarachik and is characterized by the negativeness of Lya- 
punov exponents. In general, it is not possible to determine this exponents 
explicitly. Many numerical approaches have been proposed, which generally 
used the simulation of the stochastic trajectories. In the present paper, we 
study a stochastic dynamical system that is used in economy, in describing 
a Counot duopoly game. 

In 1838, Cournot introduced the first formal theory of oligopoly, which 
treated the case of naive expectations, where each player assumes the last 
values taken by the competitors without estimation of their future reactions 
[5j. Recently, a lot of articles have shown that the Cournot model may lead 
to a cyclic or chaotic behavior [3], [1], [9], [11], [12], [TBj. Also, in [14j, Rosser 
reviews the development of the theory of complex oligopoly dynamics. 

In the present paper we have studied a stochastic Cournot economic game. 
In Section 2 we present the Lyapunov exponent and stability in stochastic 2d 
dynamical structures. Section 3 describes the Lyapunov function method for 
the stochastic stability analysis. Section 4 studies the Lyapunov exponent for 
an economic game with stochastic dynamics. The Lyapunov function method 
for the stochastic game is given in Section 5. Some numerical simulations 
are done in Section 6. Finally, Section 7 draws some conclusions. 

2 The Lyapunov exponent and stability in 
stochastic 2d dynamical structures. 

Let (O, 9^, y) be a probability space 0. It is assumed that the cj— algebra 
9^ is a filtration that is, 3" is generated by a family of cr— algebra 3^t{t > 0) 
such that 

S'sCS'tCS', Vs<t,s,tG/, 
where I =[0,T],T e {0,oo). 



2 



Let {x{t) = {xi{t),X2{t))}t>o be a stochastic process. The system of ltd 
equations: 

dxi{t, Lj) = fi{t, x{t, uj))dt + gi{x{t, uj))dw{t, u),i = 1,2, (1) 
with the initial condition x(0) = xq is written as: 

Xi{t,u) = Xio{u) + fi{x{s,uj))ds + gi{x{s,uj))dw{s,uj),i = 1,2, (2) 
Jo Jo 

for almost all u ^ Q and for each t > 0, where fi{x) is a drift function, 
gi{x) is a diffusion function, fi{x{s))ds, i = 1, 2 is a Riemann integral and 
/o giix{s))dw{s) is an ltd integral. It is assumed that and gi, 2 = 1,2 satisfy 
the conditions of existence of solution for this SDE with initial condition 
x(0) = Co e IR". 

Let Xq = {xio, X20) G IR^ be a solution of the system: 

Mxo)=0,t = l,2. (3) 

The functions gi,i = 1,2 are chosen so that: 

gi{xo) = 0,z = 1,2. 

In what follows, we consider: 

2 

9iix) =^hj{xj - Xoj),i = 1,2, 

where bij G IR, i, j = 1, 2. 

The linearized system of in xq, is given by: 



where 



X{t) = [ AX{s)ds + [ BX{s)dw{s) 
Jo Jo 



^{t)={ ^i(J'^)yA= f ],B 

' U2{t,Uj) J \ 021 ^22 ' 

..-Ml h 



Jll Ui2 
^21 ^22 



The Oseledec multiplicative ergodic theorem [TU] asserts the existence of 2 
non-random Lyapunov exponents A2 < Ai = A. The top Lyapunov exponent 

is given by: 

A = lim sup log a/ ui{ty + U2{tY- 

Applying the change to polar coordinates: 

x{t) = r{t)cose{t),y{t) = r{t)sin9{t) 
by writing the Ito formula for 

hi{ui,U2) = ]-\og{ul + ul) = \og{r),h2{ui,U2) = arctg{—) =9. 
we get: 

Proposition 1 [8]. The formulas 

log =£?l(^(^)) + ^(?4(^?(^))^-?2(^^(^))^)rf3+^\2(^(^))d^(^), (4) 

e{t)=9{0)+ [ q3{e{s))-q2{9{s)qMs))ds+ [ qMs))dw{s), (5) 
Jo Jo 

hold, where 

qi{9) = auCos'^{9) + (ai2 + a2i)cos9 sm9 + a22sin'^9, 
q^{9) = buCos'^{9) + (612 + b2i)cos9 sin9 + b22siri^9, 
qz{9) = a2iCos'^{9) + {022 — au)cos9 sm9 — ai2sin'^9, 
q4{9) = h2icos^{9) + (^22 — f'ii)cos^ sin ^ — bi2sin^9. 

As the expectation of the Ito stochastic integral is null 



(6) 



E [ q2{9{s))dw{s) = 0, 
Jo 



the Lyapunov exponent is given by: 



A = lim - lor ^ ^ 



t^oo t \r{0) 



lim [\q,{9{s)) + l{qMs))r-q2i9{s)))ds. 
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Applying the Oseledec theorem, if r{t) is ergodic, we get: 
A = l\qi{e) + ^iq,{ef - q2{9)))p{e)d9, 

where p{9) is the density of probabihty of the process 6. 

An approximation of this density is calculated by solving the Fokker- 
Planck equation. 

The Fokker-Planck (FPE) equation associated with equation ([5]) for p = 
pit, 6) is 

% + ^((^3(^?) - q2{e)qmP) - l-^2iQM'p) = 0. (7) 

From ([7]), it results that the solution p{9) of the FPE is a solution of the 
following first order equation: 

{-q3{0) + qi{0)qm + q2{0)95{9))p{0) + lqmy{0)=Po, (8) 

where p'(0) = and 
dO 

q^^e) = -(6i2 + h2i)sin2d - (622 - &ii)cos20. 

Proposition 2 |^. If qi{0) 7^ 0, the solution of the equation ^ is given 
by: 

where k is determined by the normality condition 

/ p{e)de = 1 

Jo 

and 

D{2n) - 1 
lo" D{u)du 

The function D is given by: 

Die) = exp(-2 r ^3H-g2Hg4(n) -g4Hg5H^^^ 

Jo q^iuy 
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A numerical solution of the phase distribution could be given by a simple 
backward difference scheme. 

TT 

We consider G 1R+, ^ ~ ^'^'^ 

qiii) = ail cos^{ih) + {au + 021) cos{ih) sin(i/i) + 022 siii^^ih), 
q2{i) = bii cos^(i/i) + {hi2 + 621) cos{ih) sm{ih) + 622 sin^(i/i), 
^3(2) = a2i cos^(z/i) + (022 — flu) cos(i/i) sin(2/i) — ausin'^iih), 
^4(2) = 621 cos^(i/i) + (622 — &11) cos{ih) sin(i/i) — 612 sin^(i/i), 
Q^ii) = -{bu + &21) sin(2i/i) - (622 - bu)cos{2ih),i = 0, 

The function p{i), i = 0, N is given by the following relations: 

pW = (p(o) + *(!M^)fW 

where 

2/1 

2h{-q3{i) + ^2(094(0 + Q4ii)Q5{i)) + ?4(^)' 
The Lyapunov exponent is A = A(A), where 

^ 1 

1=0 

From Proposition [2] we obtain: 
Proposition 3 // the matrix B is given by: 

bii = a, bi2 = -(3, 621 = /3, ^22 = a 

then 

k 1 1 1 

Pi^) = ^exp{-^((a2i - ai2 - af3)6 + -(an - 022) cos 26* + -(021- 

-ai2)sin26')} 



Jq"" exp{^{{a2i-ai2-al3)9 + ]-{aii-a22) cos 26*+^ (021-012) sin26')d6' 
A = ^(aii + a22 + - a^) + ^(on - 022)02 + ^(021 + 012)52, 

where 

C2= cos{2e)p{e)d9, S2= sin{2e)p{e)de. 
Jo Jo 



3 The Lyapunov function method for the stochas- 
tic stability analysis. 

We consider the system of stochastic equations, SDE, given by: 

dx,{t) = h{x{t))dt + g,{x{t))dB\t), 1 = 1,2, (9) 

where x{t) = x{t,uj) and B^{t), B^{t) are Wiener processes. 

Let V : D = {0, oo) x ^ ]R be a continuous function with respect to 
the first variable and a class function with respect to the other variables. 

Let: 

, dV(t,x) ^^^^dV(t,x) ^ ^d^V(t,x) 

i=l i=l ■' 

be a differential operator. 

We suppose that Xe = is the stationary state of ([9]), that means: 

/i(0) = (7ia(0) = 0, z = l,2,« = l,2. 

The theorem which gives us conditions for the stability of the trivial 
solution = in the terms of the Lyapunov function is: 

Theorem 4 IW^ Under the above conditions, if there is a function V : D ^ 
IR and two continuous functions u, v : 1R+ —>■ 1R+ and k > so that for 
\\x\ \ < k the relation: 

u{\\x\\) < V{t,x) < v{\\x\\) 

holds, then: 

(i) if LV{t,x) < 0, X G (Ojfc), then solution of ^ Xe = Q is stable in 
probability; 

(a) if there is a continuous function c : 1R+ — > 1R+ so that LV{t,x) < 
— c(||z||) then solution Xe = of ^ is asymptotically stable. 

Let V : D = {0, oo) x IR^ — > IR be a continuous function with respect to 
the first variable and a class function with respect to the other variables. 

The theorem that gives us the exponential p-stability condition of the 
trivial solution (fTT!) is: 
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Theorem 5 IT^ If function V satisfies the following inequalities: 

ki\\u\\P < V{t,x) < k2\\u\\'' 

LV{t,u) < -h\\u\\P, ki >0,p>0,i = 1,2, 

then the trivial solution of / flT]) is exponential p-stable for t > 0. 
For the concrete problems the following theorem is used: 

Theorem 6 /6/ If function V satisfies the following inequality: 
(i) LV{u) < 0, then the trivial solution is stable in probability; 
(a) LV{u) < —c{\\u\\), where c : 1R_|_ 1R+ is a continuous function, 

then the trivial solution is asymptotically stable; 

(Hi) LV{u) < —q^Qq, where Q is a symmetric matrix positively defined, 

then the trivial solution is mean square stable. 

In general, the functions /«, gia, i = 1, 2, a = 1, 2 are nonlinear functions 
and the above theorem is difficult to use. Therefore, the linearization method 
of system ([H]), in the neighborhood of the equilibrium point is used. 

The linearized stochastic differential system SDEL of ^ is given by: 



dui{t) = (aii-ui(t) + ai2U2{t))dt + (bii'Ui(t) + bi2U2{t))dB^{t) 
du2{t) = {a2iui{t) + a22U2{t))dt + (foai^^i W + b22U2{t))dB^{t). 

For (fTTl) expression LV is given by: 

dV dV 

LV = (aiiMi + 012^2) 7^ \- (a2iMi + 022^^2) 7^ 

OUi ou 



2 

1 d'^V d'^V 

-[{biiUi + hi2U2f-^ + (621M1 + b22U2f-^] 



(11) 



(12) 



4 The Lyapunov exponent for an economic 
game with stochastic dynamics. 

Two firms enter the market with a homogenous consumption product. 
The elements which describe the model are: the quantities which enter the 
market from the two firms > 0, z = 1,2; the inverse demand function 
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p : ]R_|_ R_|_ {p is a derivable function with p' (x) < 0, limp{x) = 0, 

X — >-ai 

limp{x) = bi, (tti G M, 61 G M); the cost functions Cj : M+ ^ M+ ( Cj are 
derivable functions with Q (a^j) > 0, > 0, i = 1, 2 ). 

In our study we consider p(x) = —,x > and CAxi) = CiXi + di,i = 1, 2. 

The mathematical model of the stochastic dynamic economic game is 
described by the stochastic system of equations: 

/■* X2(s) /■* 

xi{t) = xi{0)+ki {- — — ——-ci)ds+ {bnXi{s) + bi2X2{s) + -f^)dw{s) 

Jo {xi{s) + X2{s)y Jo 

X2{t) = X'lO)+k2 (7 — ^—^ -—--C2)ds+ {b2iXi{s) + b22X2{s) + -f2)dw{s) 

Jo {Xl{s)+X2{s)y Jo 

(13) 

where bij G H, i,j = 1,2, ki > 0,k2> 0, Xi{t) = Xi(t,u), i = 1,2. 

_ bnC2 + &12C1 _ b2lC2 + b22Ci 

(Cl + C2)2 '^2" (Ci + C2)2 • 

For bij = 0, i,j = 1,2 model (fT3|) is reduced to the classical model of the 
economic game [3], [9]. 

The system of stochastic equations (fT3!l . has the form ([2]) from section 2, 
where: 

Xa) = 7 - Ci, X2) = 611X1 + 612X2 + 7i, 

(Xi + X2)"^ 

/2(Xi, X2) = 7 ^ - C2, ^2(2^1, 2:2) = 621X1 + 6222^2 + 72- 

[Xi + X2) 

Applying the results from section 2, we have: 
Proposition 7 (i) The stationary state of (SDE) ( f73]) is given by: 

_ C2 Ci 

^10 — 7 ; ^1X20 



[Cl+C2)2' (Ci+C2)2' 

(a) The elements of the matrix A, which characterize linearized equation 
T^) in (xio,X2o) are: 

au = -2fcici(ci + C2), ai2 = -ki{cl - cl) 
021 = k2{cl - C2), 022 = -2A;2C2(ci + C2); 
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(in) The roots of the characteristic equation: 

- («11 + 0-22)1^ + CLllO'22 - «12a21 = (14) 

have the real part: 

^e(/ii_2) = -{kici + k2C2){ci + C2); 

(iv) If bii = a, bi2 = —(3, 621 = /5, ^22 = 7^ 0, t/ien t/ie Lyapunov 
coefficient of (SDE) ^ is: 

\=-{kiCi + /C2C2)(Ci + C2) + ^(/?^ - a^) - (/CiCi - k2C2){Ci + C2)£'2 + 

1 ^ (15) 

+ -(A:2-A:i)(c2-c^)^2 

/'27r r2Ti 

D2= cos{2e)p{e)d9,E2= / sin{2e)p{e)de 
Jo Jo 

and 

p{e) = kg{e),k 



g{9) = -^exp{-^((A;i + k2){cl - cl) + a(3)9 - {k^ci - k2C2){ci + C2) cos(26') + 
+ ]^{ki + k2){c\-cl)sin{2B)}. 



5 Numerical Simulations. 

We have done the numerical simulations using a program in Maple 12. 
For Cl = 0.2, C2 = 2, /ci = 0.2, k2 = 0.4, /3 = 2, in figure 1 is displayed 
(a, A(a)), where A(a) is given by f|T5l) . For a G (—00, —1.2) U (1.1, cx)), the 
Lyapunov exponent is negative, then (SDE) has an asymptotically stable 
stationary state. For a G (—1.2, 1.1), the Lyapunov exponent is positive and 
(SDE) has an asymptotically unstable stationary state. 
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Fig 1. (a,A(a)) 



/ 0.4- 
/ 0.2- 




-2 / ' ^ 
/ -0.2- 

/ -0.4- 

/ -0,6- 


l\ 2 


/ "0.8- 
/ "1.0" 
h1.2- 





If /3 is a real parameter and a = 2, in figure 2 we have: A(/3)). 

Fig 2. {p,xm 




For (3 G (— oo, — 2.6) U (2.6, oo) the Lyapunov exponent is positive and 
(SDK) has an asymptotically unstable stationary state. For /?(— 2.6,2.6) the 
Lyapunov exponent is negative and (SDE)equation has an asymptotically 
stable stationary state. 
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The Euler second order scheme for (SDE) (fT^ is given by: 

X2(n) 



Xi{n + 1) = Xi(n) + h 



+ X2{n)) 



2 - Ci ) + (611X1 (n) + bi2X2in) + 7 J 



■ G{n) + f>ii(feiixi(n) + h2X2in) + 7i)^^^^^7^ + 



2xi{n)x2{n) 
{xi{n) + X2{n))^ 



X2{n) 
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{xi{n) + X2(n))2 
2z2(n) 
(xi(n) + X2{n)y 



\ n / \ 7 / \ N Xi(n)x2(n) 
- ci + (feiiXiH + buX2in) + 7i) . . ^ \ \.J — + 



(xi(n) + X2(n))3 2 



)(6iiXi(n) + 6i2X2(n) + 71) 



hGjn) 



X2{n + 1) = X2{n) + /i 



Xi(n) 



(xi(n) + X2(n))2 

^(n)^ - h 



■ G{n) + 622(621X1(72) + h22X2{n) + 73) 



C2 ) + (6212:1 (n) + 622X2(71) + 72) 

2x1 (72)^2 (n) 
^ (xi (77.) + 0:2(77)) 3 



2:1(77) \ /, / X , /X X Xi (77)^2(77) , /7^ 

^ ^ - C2 + (621X1(77) + 622X2(77) + 72) . / x \\xJ — + 



(621- 



(Xi(77) + X2(r7))2 
2X1(77) 



(Xi(77) + X2 ( 77))3 2 



(Xi(r7) + X2(77))3 



) (621X1 (77) + 622X2(77) + 72) 



/7G(77) 



where G'(7t,) = 77;((77 + l)h) — w{nh), n = 1,2,..., and Xj(77) = Xi{nh,uj), 
7 = 1,2. 

In figures 3 and 4 the orbits: (77, Xi(77, a;)) for (SDE) and (77, Xi(77)) for 
(ODE) are displayed: 



Fig 3. (77, xi(77, uj)) 



Fig 4. (77, Xi(77)) 



0.45- 
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0.41- 
0.40-^ 
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200 400 600 800 1000 
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In figures 5 and 6 the orbits: {n,X2{n,uj)) for (SDE) and {n,X2{n)) for 
(ODE) are displayed: 



Fig 5. {n,X2{n,(jj)) 



Fig 6. (n, X2{n)) 



20U \ 400 600 SOO 1000 



200 400 600 800 1000 



In figures 7 and 8 the orbits: (xi(n, u), X2{n, u)) for (SDE) and {xi{n),X2{n)) 
for (ODE) are displayed: 



Fig 7. {xi{n,u),X2{n,uj)) Fig 8. {xi{n) , X2{n)) 



0.36 0.37 0.3^ 0.3? 0.40 :0.41 0.42 0,43 0.44-0.45 



0,400 0,402 0,404 0,406 0,40S 0.410 0,412 



6 The Lyapunov function method for the stochas- 
tic economic game. 

Theorem [6] is used for the analysis of the stability with the help of the 
Lyapunov function. 
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Let V : D = {[0,(yo) X H^} ^ H be the function given by: 

where Ui > 0,i = 1, 2. 

Using formula ( fT2l) for the hnearized system of ( fT3l) : 

dui{t) = (aii-ui(t) + ai2U2{t))dt + {hiiUi{t) + hi2U2{t))dB{t) 
du2{t) = {a2iUi{t) + a22U2{t))dt + {b2iUi{t) + h22U2{t))dB{t), 

we obtain: 

Ly(-U(t)) =(011111 + 012^2)0^1^1 + (a2l'Ui + a22tt2)w2'W2 + 

^[(fellMl + 612^2)^^1 + (&21«1 + &22M2)^tU2] = 

11 11 

= (aiiu;i + -bliUJi + -hl^uj2)ul + (022^^2 + 2^12^1 + 2^22^2)^1+ 



(16) 



+ (ai2Wi + a2iu;2 + 6iibi2u;i + 62i622u;2)'Ui'U2- 



(17) 



If 



a2i + 621&22 , u u / A 

Ai = ■ , , ai2 + 611612 7^ O.cui = -Aiu;2, 

ai2 + 611612 

11 11 
qi = {an + ds-bu)Ai - -6^^, ^2 = -^22 - ^^22 + 2^12^1' 

then from (ITTI) and (fT8|) we get: 

LV{u) = -qiUJ2u\ - q'2W2'U2- 

Form the above relations and Theorem [6] we obtain: 
Proposition 8 Ifbij,i,j = 1,2 satisfy the relations: 

ai2 + 611612 ^0,Ai < 0, qi >0,q2> 0, 
i/ien the trivial solution of [TS\) is mean square stable. 



(18) 
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7 Conclusions. 



In the present paper we investigate an economic game with stochastic 
dynamics. We focus on a particular game and determine the Lyapunov expo- 
nent for the stochastic system of equations that describes the mathematical 
model and the Lyapunov function for the analysis of the mean square sta- 
bility. The calculation of the top Lyapunov exponent allows us to decide 
whether a stochastic system is stable or not. Using a program in Maple 12, 
we display the Lyapunov exponent and the system orbits. Conditions for the 
solution of the stochastic game to be asymptotically mean square stable are 
estabhshed. 
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